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. , Abstract 

'NT 

o 

O ' We study the phase diagram, ground state properties, and excitation gaps of a 



frustrated spin-Peierls chain, i.e. a spin-Peierls chain with both nearest neigh- 
bor exchange Ji and next nearest neighbor exchange J2. The phase diagram 



' is calculated using a bosonization Renormalization-Group (KG) method. We 



propose a (bosonic) bond-operator mean field approximation to calculate the 
ground state properties and excitation gaps. 

I. INTRODUCTION 

A theoretical understanding of the frustrated spin-Peierls chains {i.e. frustrated Ji — J2 
Heisenberg chains with spin-lattice interaction) is important for a number of reasons: First, 
there are a number of quasi-one- dimensional systems which can be described by a spin-Peierls 
model at low temperatures, and in some of these systems, such as CuGeO^ materialJll and 
some MX material^!, the competition between nearest neighbor Ji and next nearest neighbor 
J2 plays an important role. Second, the frustrated spin-Peierls model poses interesting 
theoretical problems. The frustrated Ji — J2 Heisenberg model has been studied in the 
literatureiHB. MajumdaJl showed that at J2/ Ji = 1/2 the ground state is dimerized and the 
exact ground state wave function is 
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^i = [l,2][3,4]---[2iV-l,2iV] 

^2 = [2Ar,l][2,3]---[2iV-2,2Ar-l], (1) 

where denotes the normahzed singlet combination of the spins on sites i and j. Later 
it was shown by HaldaneS that dimerization exists at J2/J1 > Jc ~ 1/6- For J2/ Ji < Jci 
the system is in a gapless spin hquid phase. This resuh was confirmed by numerical studies 
on finite chainsBi. The critical value Jc found by numerical studiesi'i is Jc ~ 0.3. With 
spin-lattice interaction, there is dimerization even without J2 frustration due to the Peierls 
instability. The interplay of J2 frustration and spin-Peierls mechanism is an interesting 
subject. 

In this work, we consider the frustrated spin-Peierls system described by the Hamiltonian: 
H = Y.^JiA ■ S^+l + J2S, ■ S,+2) + \ Y.^Miu] + Kul), (2) 

where Ji^i = Ji — PuUi- In Sec. II, we recalculate the critical value Jc for the Ji — J2 
Heisenberg model using Haldane's method and obtain Jc ~ 0.4. We also study the RG 
equations for the spin-Peierls system in the adiabatic limit and find that the spin-lattice 
interaction is always relevant and opens a finite gap. In Sec. Ill, we use a bond-operator 
mean field approximation to calculate the dimerization magnitude, the ground state energy 
and the excitation gap. We find that this approximation can produce a very accurate ground 
state energy for the Ji — J2 model for J2/J1 < 0.6, and the ground state energy and gap 
calculated by this solution is exact at J2/ Ji = 1/2. Since our approximation is based on the 
dimerization of the system, we belive our mean-field approximation should be reliable for 
the spin-Peierls system. Our results are summarized in Sec. IV. 



II. SCALING THEORY 



A. Frustrated Ji — J2 Model 

Here we first discuss the frustrated Ji — J2 model studied by Haldanei for the S = ^ 
spin system 

Hs = Ji{SiSf^i + SfSf^i + 'yS-S-j^^i) + J2Si ■ Si+2, (3) 

where the anisotropy parameter 7 is introduced. 

Using the Jordan- Wigner transformation, Hamiltonian @ transforms into the spinless 
Fermion Hamiltonian 
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Hs = Yl -y (4ci+i + H.C.) + Jaim - l/2)(n,+i - 1/2) + 
Mn, - l/2)(n,+2 - 1/2) + [cl(n,+i - 1/2)q+2 + i^-C] 



(4) 



In the continuum limit, this Hamiltonian becomes 



Hj{JiaQ) = J dx{ii)\Vi)+ - itplV^-) + g^ip+p^ + p-p-) + 2g2{p+p-) 

-g3{ip+ix)ip_{x)ijjl{x + ao)tp_{x + ao) + H.C), 



(5) 



where g^ = ■y + 2J2/J1, (72 = 27 + 2J2/J1, and gs = ■y — 3J2/J1. Here our coefficients gi are 
different from Haldane'si. The crucial difference is in the g3 term. The difference comes from 
the term ip\_{x)tp_{x + 2ao)ilj\^{x + ao)ilj^{x + a^) . This term is written as {tp\.'Vip\.)ip-'Vtp- in 
Refl, which has opposite sign to the term from il!\_{x)ip^{x)il!\_{x+2aQ)ip-{x+2aQ) . However, 
by direct bosonizationi0, both expresses contribute negative Umklapp term: cos(40). Here 
is the boson phase field: 



where N± are the number operators. Thus we have g^ = 7— 3 J2/ Ji instead of g^ = 7— 6 J2/ Ji- 



7 ~ 1 are identified by HaldaneS as Neel state, gapless spin liquid, and dimer state (with 
a finite gap). The scaling digram shown in Fig.(|l]) is calculated numerically by integrating 
the scaling equations. The critical value of J2/J1 is calculated numerically by locating the 
intersection of the phase boundary curve and curve l/Kp = [(7r+ 13/37 — 4/3(y'3)/(7r — 7)]^/^. 
The critical value we obtained is Jc ~ 0.447 instead of 0.16. The critical value calculated 
by Density-Matrix-Renormalization-Group method using 300 sitesi is Jc = 0.298. Since we 
only used the lowest order scaling equations and the critical point of J2/J1 is not close to 
1/Kp = 2.0, the value of Jc calculated here is not accurate. The true critical value of Jc 
should be close to the numerical valud^i 0.3. 




[P+{P) + P-ip)] - (N^ + N-) 



(6) 
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FIGURES 




FIG. 1. Scaling Diagram for Ji — J2 Model 



B. Frustrated Spin-Peierls Model 

The lattice part of the Hamiltonian (|^) can be written as follows (after rescaling u and 
going to the continuum limit): 



2 'luU2kf 



cos(20)), 



(7) 



where a is the cutoff and is the boson field defined in Eq.(^. We have use the following 
equation to derive Eq.(0): 



with the definition 



(8) 



0±(a;) = - ((p{x) =F J do(f){x)dx'j . 
In terms of the boson field (j){x) the total Hamiltonian is expressed by 



(9) 



H = Hq + Hp + Hj, 



(10) 



where 



Hr 



1 

2^ 



dx 



\jdx 
(27ra)2 



1 

Pi 



and 



1 

dx cos(40) H 

Tca 



TT 



dx cos{2(j))u2kf 



(11) 



+ 9a- 92 

+ 9i + 92 



TT 
1/2 



1/2 



(12) 



With the Fermion-phonon interaction, the scahng theory is much more comphcated. The 
analysis is only possible in high phonon frequency limit A <^ uj2kf or in the adiabatic limit 
A ^ uj2kf, where A is the gap in the spinless Fermion spectrum. In this paper, we will 
concentrate on the adiabatic limit, which is valid for many compounds. For completeness, 
we first discuss the high phonon frequency limit. 

(a) uj2kf 00 

In this limit, the phonon degrees of freedom can be integrated out completely 
and the result is to introduce an effective attractive electron-electron interaction: 
— 72//i:/(27ra)2cos(20). Then the phonon effect is to renormalize the Umklapp interaction: 



93 = 93 - I'^/k. 



(13) 



The phase boundary between the spin liquid and dimerized state is changed accordingly. 
For weak electron-phonon coupling <^ 1, the dimerization critical point changes from 
Jc ^ Jc — For strong electron-phonon coupling 'j'^/k > 1, the system is always 

dimerized. 

(b) finite high phonon frequency limit uj2kf ^ A 

For high phonon frequency u!2kf ^ A and weak electron-phonon interactions, we can still 
integrate out the phonon degree of freedom and obtain a phonon-induced retarded attractive 
electron-electron interaction0lli: 



S93 



cos(40) = -7V^ dt dt'ip{x,t)^/j{x,t)Do{t - t')^/j{x,t')^p{x,t') 



iVk 
{2-KaY 



cos(40) - (7Vfc) 



n=l 



Qn 

dt^ 



J{x, t) 



2n 



(14) 



where 
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J{x,t) = ip{x,t)ip{x,t) 

= ilj\.{x,t)i/j-{x,t) + ilj'L{x,t)ilj+{x,t) 

and -Do(t) = uj /2 exp{—i\t\uj) is the bare phonon Green's function. All the terms containing 
explicit derivatives are superficially irrelevant operators. So for finite high frequency phonons 
and weak electron-phonon interaction, the retardation effects will be scaled out at /* ~ 
Vf/uj2kf or Ef{l*) ~ co'2fcj:& The critical theory is governed by the limit uo2kf — * oo. 
(c) adiabatic limit U2kf — ^ 

In the adiabatic limit, we should "integrate" out the electronic degree of freedom first 
and then study the effective phonon system. To this end, we can use the following strategy to 
study the spin-Peierls system in the mean- field approximation for the phonon effect: we can 
first study the electronic Hamiltonian by assuming a mean lattice dimerization U2kf = Uq] 
then we can calculate uo by minimizing the electron-lattice total energy. In this adiabatic 
limit, the electronic Hamiltonian we want to study is: 



H = — I dx 

2tt 



+ / dxcoM) + ^ / cixcos(20). (15) 

[2'KaY J Tia J 

In lowest order, the scaling equation can be derived using the Coulomb gas analog y00. 



The resulting equations are: 



dY. 1 
d\na ■ 2 



' -2Y,{l-2Kp) + -Yl (16) 



,,^,-yA2-Kp^y,) (17) 

"^^^ - ^iYl + Y^% (18) 



where Yg = g^/n and Yph = 'juUoa. The scaling diagram Fig.(|D following from Eq.(|l6D-(|l8D 
is calculated numerically by standard ordinary differential equations integration. 
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FIG. 2. Scaling Diagram 



In the Coulomb gas analogy, the system is equivalent to a Coulomb gas of four kinds 
of particles {i.e. two fugacities Yg and Yph) with charges (±eo,±2eo). Here /3eo = {26)~^ . 
The activities Aj of these particles are exp(— /?Ai) = Yg/{2aY for (±2eo) particles and 
exp(— /9A2) = Yph/ {2a) for (±eo) particles. The transition from the gapless phase to the 
gapped phase is of Kosterlitz -ThoulesS type. If there is no coupling between these two 
fugacities, the two coupling constants Yph and Yg define two separate transition "temper- 
atures" approximately at \Yg\ = |1 — and \Yph\ = |1 — 1/(46*) |, corresponding to the 
breaking of pairs (— 2eo, 2eo) and (— eo, cq), respectively. For each of the fugacities, the scal- 
ing trajectories are well known0'0, see Fig. |l| for the scaling diagram of 1/Kp — Yg, the 
scaling diagram of 1/Kp — Yph can be obtained by rescaling 1/4 of 1/Kp axis in Fig. |1|. 
The coupling between the fagucites contributes to the scaling equations in two respects: (1) 
Kp is renormahzed by Yph and Yg simultaneously; see Eq.(|T^); (2) the close "pairs" with 
nonzero net charge will contribute to the renormalization of Yph and Yg. Physically, it is 
easy to see that the transition from gapless to gapped phases is controlled by the unbinding 
of (eo, —Co) pairs, since the unbinding of (cq, — Cq) pairs has a lower transition temperature. 
So the scaling process is dominated by the Yph scaling (for initial Yph 7^ 0) (see Fig. H): the 
separatrix line is through the point {Yph = 0,Yg = 0, l/Kp = 1/2). The gapless spin liquid 
phase exists only at 1/Kp < 1/2. Since for general value of Ji and J2, 1/Kp > 1, the system 
is always in the gapped phase. 

Since we have used the bosonization procedure, Qi and 7uMo are treated as small param- 
eters. Consequently, we can not expect that the above scaling theory is correct for large 
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or large J2I J\- Most importantly, since we used a continuum model, "lattice-effects" can 
not be recovered. From analytical results! we know that in the special limit J2/ Ji = 1/2, 
the ground state is a perfect dimer configuration which should be true only for a discrete 
model. Intuitively, we may expect that if we keep J2 finite and Ji 0, the system behaves 
as two decoupled gapless spin liquids. This phase does not occur in the above scahng phase 
diagram. Nevertheless, from the scaling theory, we can know that for the Ji — J2 Heisenberg 
model, there is a gapped phase at J2/J1 > Jc (but not J2/ Ji — ^ C)o), and from the exact re- 
sults at J2/ Ji = 1/2, the ground state in this gapped phase is dimerized. For the frustrated 
spin-Peierls model (in the adiabatic limit), from the scaling theory, there is always a gap 
due to spin lattice coupling. 



III. BOND-OPERATOR MEAN FIELD SOLUTION 

From Ref.i, we know the ground states of the J1 — J2 model at J2/ Ji = 1/2 are degenerate 
dimer states with wave functions 

^i = [l,2][3,4]---[2iV-l,2iV] 

^2 = [2iV,l][2,3]---[2A^-2,2A^-l]. (19) 

The continuum limit discussed above showed that the gapped state will exist at J2/ Ji > Jc ^ 
0.3. Thus we can expect these gapped states to be dimer states similar to the J2/J1 = 1/2 
case. Since we know that in this region the low energy states tend to have dimer valence 
bonds, we will use the bond-operator representation of the spins first proposed by Sachdev 
and Bhatt for the 2D mo del§. 

1 2 3 4 2i-1 2i 2N-1 2N 



FIG. 3. Dimerization of the Chain 

The bond-operator representation depends on the dimerization of the valence bond. Let 
us choose a dimer pattern as in Fig ^ For each pair of spins connected in the dimer pattern, 
we can use four bond spin operators which are defined by the singlet |s) and triplet |tig{i,2,3}) 
states: 

k) = ^^|o) = ^(in)-|iT)) 
|ti) = tl|o) = ^(|TT)-Ui)) 
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N = 4io) = -^(in) + ui)) 

l^3) = 4|o) = ^(in) + UT)) (20) 

The the representations of the four bond operators s„ and tn,ie{i,2,3} is a faithful repre- 
sentation of the original spin operators and these bond operators can be treated as boson 

— * 

operators. Their relation to the spin operators Sn is as follows: 



^2n ~ 9 ( ^n^'>i,a in,a^n i^afii 



2 

Using the bond operators, the Hamiltonian becomes: 

H^Ho + Hi + H2 + Hs, (22) 

where 

3 1 - - 1 

Ho = l](^l + u){--slsn + -riijn,a) + -NKu^ - Y.^&n + tl,a4,a - 1) 



4 

Hl=J2 + " + 2^2)(tl,a4+l,a^nSn+l + ^l,a4+l,a)s|,+iS„ + H.C.) 

n ^ 

-f^2 = X] 7("^1 - 1i + 2J2)(-4,a^+l,aWl,/34,/3 + ^l,Q:C+l,/34+l,a4,/3) 



n 

H3 — -[Ji — U)eal3'y{itjiatn+l,l3tn+l,'ySn + ^^n+l,a!^n,^^n,7^n+l + H.C). (23) 

Here wc have used a mean field approximation for the phonon part. Because of the 
— (3/4)sJjS„ term in Hq, it is clear that the s bosons will condense: = Sq. The terms 
in Hi and H2 suggest a nonzero value of (iaia)- As J2 increases, the ti bosons will also 
condense. Whether t, condense or not, we can omit H3 in the mean field approximation. 
Using this procedure, the Hamiltonian becomes: 

H = N[-^{Ji + u)sl -fisl + fi + Ku^] + ^-{%h + Ji - n)B'' - ^-{23^ + Ji - u)J^ 

+(J2 + Ji/2 - K/2)AcosA;]iXJfe,„ + ^[7(2^2 ^u- Ji)B cos k]t\Jl,^^ + H.C. 

k ^ 

+ Y^[j{Ji + u) - II + {J2 + u/2 - Ji/2)slcosk. (24) 

k ^ 

The self-consistent equations (at zero temperature) are: 

A — sq -\ — ^ cos k 
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R ^ 3 ^ A, 



(25) 



and 



efc = 


{Ji + u)/2 + ^{2J2 + u- Ji] 
-(2J2 — u + Ji)A cos 


Sq cos k 


-i(2J2 + M- Ji)(A + 5) + 


Afe = 


^[{2J2 + Ji-u)B-{J2 + u 


- Ji)sl] 


cos A; 


2Ku = 




1 

4 




Ek = 


^4 - Ai. 







(26) 



The self-consistent equations ( pSf ) and (^) can be solved iteratively. The numerical 
results are summarized in Figi 
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FIG. 4. Ground state energy for J1 — J2 Heisenberg model (solid line) and frustrated spin-Peierls 
model (dashed line) at K = 100. 
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FIG. 5. Comparison of ground state energy from our mean- field approximation (solid line) and 
that from DMRG calculation (diamond) for Ji — J2 Heisenberg model. 

In the limit uq = 0, the above Hamiltonian describes the Ji — J2 model studied 
earlieri'ifl'i. We consider this case first. From equations (p^-(|26D, it is easy to show that at 
J2/ ^1 = 1/2 the ground is at sq = 1 and A = B = 0, i.e. the ground state is perfectly dimer- 
ized. The ground state energy is — 0.375Ji and the gap value is 0.25Ji, which are exactSi. 
The ground state energy we calculated (Fig. ^) agrees with exact numerical resultsB'i for 
J2/ t/i ~ — 0.59 Ji (see Fig. | for comparison to the calculation in Ref.i). This is surprising 
since we do not expect dimerization for small J2 without the spin-lattice interaction. Our 
interpretation is that at small J2, although there is no dimerization, the ground state is still 
well described by fluctuating singlet valence bonds, so our mean-field bond operator approx- 
imation can give a rather accurate ground state energy. Since the spin-lattice coupling will 
enhance the dimerization, we can expect that the mean field ground state energy for the 
frustrated spin-Peierls model will also be accurate. 

The lattice dimerization magnitude Uq for the frustrated spin-Peierls model is shown in 
Fig. |. Since the lattice dimerization magnitude Uq only depends on the ground state of the 
system, we believe that this mean field result should be accurate. From Figs, ^^we can see 
that the dimerization of the lattice is enhanced by next nearest neighbor antiferromagnetic 
coupling J2. 
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FIG. 6. Lattice dimerization magnitude u for the frustrated spin-Peierls model at K = 100. 



The gap calculated here (Fig. 
for the Ji — J 2 Heisenberg model 



is also reasonable comparing to the previous results 
. However there are two features of the mean-field 
gap for the Ji — J2 Heisenberg model which are incorrect: (i) We obtained a finite gap 
even for small J2/J1 values; and (2) The gap value starts to decrease immediately after 
J2/J1 = 1/2 in contrast to the DMRG results! that the gap starts to decrease around 
J2/J1 ~ 0.7. The first incorrect feature is because we artificially assumed dimerization in the 
bond-operator representation even for small J2/Ji- For J2/J1 > 1/2, bonding between next 
nearest neighbors becomes important and the description of this kind of bonding via nearest 
neighbor bonding will not be sufficient at the mean-field level. For J2/J1 > 0.59, the ti bosons 
start to condense, and the gap decreases to zero. Although the ground state energy and 
dimerization due to spin-lattice interaction can be calculated after the t boson condensation, 
since we can not expect dimerization at large J2/J1 when next nearest neighbor valence bond 
becomes strong, the bond-operator representation is questionable in this J2/J1 region. The 
bond-operator mean-field approximation should work better for the frustrated spin-Peierls 
model, since even for J2/J1 = there is then dimerization due to Peierls instability. Our 
bond-operator approximation calculation of the gap should be reliable for J2/J1 G [0,0.59]. 
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FIG. 7. Energy gap for the for Ji — J2 Heisenberg model (solid line) and frustrated spin-Peierls 
model (dashed line) at K = 100. 



IV. SUMMARY AND ACKNOWLEDGEMENT 



In this work we have discussed the scaling behavior of a frustrated spin-Peierls system. 
The ground sate energy, the excitation gap, and lattice dimerization magnitude have been 
calculated using a bond-operator mean-field approximation. We found that the dimerization 
of the spin-Peierls chain is substantially enhanced by the next nearest neighbor antiferro- 
magnetic coupling J2. In the most interesting region, J2/J1 < 0.59, our bond-operator mean 
field solution can describe the model rather accurately. 
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